Spin correlation functions and susceptibilities in the easy-plane XXZ chain 
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We present a Green's-function theory of magnetic short-range order in the 5=1/2 easy-plane 
XXZ chain based on the projection method for the dynamic spin susceptibihty and a decoupling 
of three-spin operator products introducing vertex parameters. The longitudinal and transverse 
static susceptibilities and two-point correlation functions of arbitrary range are calculated self- 
consistently for all wavenumbers, temperatures, and anisotropy parameters —1 < A < 1. In the 
easy-plane ferromagnetic region (A < 0), the longitudinal correlators of spins at distance n change 
sign at a finite temperature To{n, A), in reasonable agreement with recent data obtained by finite- 
chain diagonalizations. The temperature dependence of the uniform static susceptibilities exhibits 
a maximum which is explained as an effect of magnetic short-range order which decreases with 
increasing temperature. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.-s 



The study of low-dimensional spin systems constitutes 
a growing field of research which is mainly due to the 
availabihty of new materials, such as the quasi-one- 
dimensional (ID) cuprates [1]. Recently, an unexpected 
quantum-classical crossover of the longitudinal spin cor- 
relation functions in the ID XXZ model 
J V^, 



n = -'^{stsj + ^stsi 



(1) 



((«,j) denote nearest-neighbor (NN) sites along the 
chain; throughout we set J — 1) was found for —1 < 
A < by means of exact diagonalization (ED) of systems 
up to 18 spins [2] and by the quantum transfer matrix 
formalism [3]. 

Motivated by those findings, in this paper we examine 
the spin correlation functions in the easy-plane region 

— 1 < A < 1 of the model (1) by an analytical approach 
based on a Green's-function projection method. This 
theory was found to provide a good description of anti- 
ferromagnetic (AFM) short-range order (SRO) in the 2D 
spatially isotropic [4] and anisotropic Heisenberg mod- 
els [5]. Moreover, for the first time, we calculate the 
full wavenumber, temperature and A dependences of the 
static transverse and longitudinal spin susceptibilities. 

To determine the dynamic spin sus- 
ceptibilities x+"('7,t^) = -{{S;^;SZq))<^ and x^^(9,w) = 

— {{Sg; SZg))uj (— TT < g < tt), defined in terms of two- 
time retarded commutator Green's functions, by the pro- 
jection method, we choose the two-operator basis Ai = 
(S'+,i5'+)-^ and A2 — {S^,iS^)'^, respectively, and con- 
sider the matrix Green's function {{A; A^))^, neglecting 
the self-energy part [4], ((A; At))^ = [cu - m'm'^]-^Tl 
with 971 = ([A, At]) and M' = {[iA, At]). We get 



M. 



(2) 



where the first spectral moments Mq 
exact expressions 



,.(1) 



are given by the 



M+-(i) = -2[C+"(1 - Acosg) + 2Cf (A - cosg)] , (3) 



iv±q 



-2Ci+"(l -cosg). 



The two-spin correlation functions C+ = {S^S^ 
C^^ — {SqS^) are calculated from 



-E 



u{l) 



2wj; 



-[l + 2p«)]e^' 



(4) 



and 



(5) 



where the Bose function p{LOq) — (e"?/^ — 1)"-^ appears 
due to the use of commutator Green's functions. The NN 
correlation functions are directly related to the internal 
energy per site e = Ci -I- ACf ^ . 

To obtain the spectra lo^ in the approximations —S^ = 

{LO+-YS+ and -S^q = (wf we take the site- 

representation and decouple the products of three spin 
operators in —Sf and — S'f along the NN sequence (i, j, I) 
introducing vertex parameters in the spirit of the scheme 
proposed by Shimahara and Takada [6] , 



S^SjSi — a2^{Sf Si )Sj 



(6) 
(7) 
(8) 



Here, and ctj appearing in ujq are attached to correla- 
tion functions between nearest and further-distant neigh- 
bors functions, respectively. We obtain 

(c^+-)2 = i[(l + 2a+-C+-)(l - Acosq) 



-A(l 



Aa+-Cni^^cosq) 



1 



+2af (A cos 25 - cos (?) 
+4a5^-Cf (cos2g- Acosg)], (9) 

«)2 = (l-cosg)[l + 2arC2+- 

-2Aaf Ci+- (1 + 2 cos g)] . (10) 

In the easy-plane region — 1 < A < 1, the theory has 

eight quantities to be determined self-consistently {Ci, 
6*2, cti, Qfj) and six self-consistency equations (5) in- 
cluding the sum rules Cq~ = 1/2 and Cq^ = 1/4. To 
get the two remaining conditions for determining the 
free a-parameters a2 we may adopt different phenomeno- 
logical choices. Following the approach by Kondo and 
Yamajii [7] for the isotropic Heisenberg chain (A = 1; 
a'j = af^), let us first consider the simple conditions 
Q!2 = a'[. In this case we do not obtain quantitatively 
satisfactory results (cf. Fig. 1); therefore, we improve 
the theory by another choice with 0:2 ^ cti- For it we 
need additional conditions. At T = 0, it is natural to 
adjust a2 to taken from the exact expressions for 
the ground-state energy £(A) [8] and the NN correlator 
Cf^ [9] (see Fig. 1). Moreover, to formulate conditions 
also at finite temperatures, we follow the reasonings of 
Ref. [6] and [4] and conjecture that both "vertex correc- 
tions" tti (T) — 1 and (^) ~" 1 have similar temperature 
dependences and vanish in the high- temperature limit. 
Correspondingly, as the simplest interpolation between 
high temperatures and T = we assume the ratio R of 
two vertex corrections as temperature independent and 
given by the ground-state value. To be specific, in calcu- 
lating x^^(<i', w;T) for all A, we assume 



and 



ar{T) - 1 



(11) 



For (q^^'jT) the additional condition analogous to 
Eq. (11) (substitution of a^^ by a^~) yields a solution 
of the self-consistency equations only in the case A > 0. 
For A < it turns out that a solution can be obtained 
using a slightly modified condition (substitution of a^^ 
in Eq. (11) by ")• That is, for x+~(g,a;;T) we require 



-(T)-l 



a 



■(T) 



Rt- 



= R+- 
af%T) - 1 < 



(12) 
(13) 



The set of self-consistency equations is solved numerically 
by Broyden's method with relative error less than 10~^. 

To discuss the rotational invariance of the theory at 
A = ±1, it is useful to perform the unitary transforma- 
tion which rotates the spins on all odd sites 171 — 21 + 1 
around the 2;-axis by the angle tt [10]. Taking the uni- 
tary operator U = Hm^'^m which transforms S„ as 
S„ = U+SnU, we get 

5^''' = e-"5^'% S^ = 5^, (14) 



(15) 



{hi) 



Since {A)^ = 
the relation 



(A)^ for any operator A, by (14) we obtain 



(9,'^) =xi (fc,w); k = q 



TT 



(16) 



For the correlation func- 
and C,^f„ = C^V,. 



and Xw(9,w) = 

tions we have . ...... ^^^.^ 

Evidently, for A = 1 the rotational symmetry is pre- 
served, i.e., C+~ = 2C^^, af~ = af^, a;+~ = w^^, and 
X"' (q,i^) = 2x^^{q,i^)- At A = —1, the rotational sym- 
metry occurs in the transformed model (15), where the 
relations X'^~{q,^) = XH~il + ^'^) = 2x^''(9,w) and 
2C^^^ = can be easily verified from Eqs. (3), (4), 
(9), and (10). 

Corresponding to the unitary equivalence of Ti. and 
Ti we denote the (paramagnetic) easy-plane region with 
— 1 < A < as ferromagnetic (FM) region (cf. Eq. (15)) 
and the region with < A < 1 as AFM region (cf. 
Eq. (1)). 

Concerning the question of magnetic long-range order 
(LRO), in the easy-plane XXZ chain there is no LRO 
which is correctly reproduced by our theory. In Ref. [5] 
the description of LRO within our scheme (mode conden- 
sation) is outlined for the spatially anisotropic Heisen- 
berg model, and the absence of LRO in the ID limit is 
shown. Considering the model (1) at ]A] < 1 and allow- 
ing for a possible finite condensation part C in the spin 
correlators [5] , a solution of our self-consistency equations 
only exists for C = 0, i.e., there is no LRO. 

In Fig. 1 the zero-temperature correlators appearing in 
the spectra (9) and (10) are plotted as functions of A. 
The NN correlators , in particular Cf^, calculated by 
the simplified theory = ai) strongly deviate from 
the exact values [8,9], especially in the vicinity of A = 0. 
The same refers to x^^(O)- Therefore, in the following 
we only show and discuss the results obtained by the im- 
proved theory [cf. Eqs. (11) to (13).] At A = 1 the rota- 
tional symmetry is visible. At the quantum critical point 
A = — 1 there occurs FM LRO [10], and a non- analytical 
behavior is observed, e.g., Um/\,^_i+ dCf^'/dA = 00 and 
limA^-i- 9Cf /9A = [8]. Moreover, at A = -1, we 
have C^~- = 1/6. This limiting behavior, however, is 

hard to obtain numerically because of the infinite slope 
of C;; as A ^ -1+. 

Figure 2 displays the static susceptibilities x''(9 = 0) "■) 
at T = vs. A. First let us emphasize the excellent 
agreement of our result for the longitudinal uniform sus- 
ceptibility with the exact result [11] over the whole A 
region (cf. inset), where X^^(O) varies by two orders of 
magnitude. That means, the longitudinal spin correla- 
tions of arbitrary range are well described by our theory. 
In the FM region, the uniform longitudinal susceptibility 
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diverges in the limit A ^ — 1 (cf. Fig. 5 (a)) indicating 
the instabihty of the paramagnetic phase against the FM 
LRO phase at A = —1 [12]. Tiie same is true for the stag- 
gered transverse susceptibihty x^Ti'^) i^^- ^ig. 5 (b)) 
which, by Eq. (16), is equivalent to the uniform trans- 
verse susceptibihty X^~(0)- 

The wavenumber dependences of the zero-temperature 
static susceptibilities x^io.) £^rc depicted in Fig. 3. For 
sufficiently low values of A in the FM region, x'''^'{q) 
shows a maximum at g = being indicative of the FM 
instability (cf. Fig. 2). Accordingly Xn~ il) ^ max- 
imum at (J = TT (Fig. 3 (b)) corresponding, by Eq. (16), 
to the maximum of X^~(^) at = 0. In the AFM re- 
gion, the longitudinal and transverse susceptibilities re- 
veal peaks at the AFM wavenumber being indicative of 
the AFM LRO at A > 1. 

Let us now discuss the finite-temperature behavior of 
the longitudinal spin correlation functions. Table I shows 
that in the AFM region we get the expected alternating 
signs of Cfj ' being characteristic of AFM SRO. By con- 
trast, in the FM region, C^^ < holds Vn below some 
characteristic temperature 7b (A) [2]. In view of the ap- 
proximations made in our theory, the analytical results 
obtained for C^^ are in reasonable agreement with the 
exact data for the 18-site chain. At fixed separation n 
and with increasing temperature or above To (A) with 
increasing separation, C^^ changes sign from negative to 
positive values. This property can be interpreted as a 
quantum-classical crossover [2] because with increasing 
temperature the system behaves more classically, i.e., it 
becomes dominated by the potential energy (longitudi- 
nal part of the Hamiltonian). It is worth emphasizing 
that the so-called "sign changing effect" in the longi- 
tudinal spin correlations found numerically [2] is repro- 
duced for the first time by an analytical theory. This 
is demonstrated in Table II. The temperatures To(n, A) 
where C^^(ro(n, A), A) = are compared with the ED 
results showing an increasing agreement with decreasing 
anisotropy parameter A. As might be expected, for the 
2D XXZ model our results for To(r, A) much better agree 
with the exact data (cf. Table 11). A detailed study of 
the 2D case is presented in Ref. [13]. 

Figure 4 shows the temperature dependence of the 
transverse correlation functions in the FM region. The 
oscillating behavior of C+~(T) (sign (—1)") and the de- 
crease of |C+~| with increasing separation n (at fixed 
temperature) and increasing temperature (at fixed n) in- 
dicates a transverse AFM SRO (but no LRO) related to 
the maximum of Xn i'l) at g = tt (cf. Fig. 3 (b)). How- 
ever, due to the equality C+J^ = (-1)"C+:^, for A < 
we prefer to denote this SRO more physically as trans- 
verse FM SRO (in the model H), where C+:^(T) can 
be read off immediately from Fig. 4. Accordingly, for 
A > 0, we get transverse AFM SRO (similar behavior of 
C+-(r) as shown in Fig. 4). 

In Fig. 5 various susceptibilities x'^{q) at g = 0, tt are 
plotted as functions of temperature. At T = 0, the di- 



vergences of xlf (0) and Xn (^) ~ (^) ^^'^ hmit 
A — > — 1, discussed above (Fig. 2), are clearly seen. In 
the region —0.3 < A < the uniform longitudinal sus- 
ceptibility reveals a maximum at Tmax which may be 
explained as follows. For T < ro(A, n) the longitudinal 
SRO, characterized by C,^^ < (cf. Table I), leads to a 
spin stiffness against the orientation of the spins along a 
homogeneous external field in .z-direction so that x^^(O) 
is suppressed. With increasing temperature the corre- 
lations become increasingly ferromagnetic, which results 
in an increase of x^'^(0;T) up to T^ax determined by 
To(n, A). With decreasing A (-1 < A < -0.3) the max- 
imum disappears, since the spin correlations are predom- 
inantly ferromagnetic. Thus, the maximum in x^^(0;T) 
for — 0.3 < A < may be understood as a combined SRO 
and sign changing effect. In the AFM region, the max- 
imum in the longitudinal susceptibility (see Fig. 5 (a)) 
obtained at all A, where T^ax shifts to higher values with 
increasing A, is due to the decrease of AFM SRO with 
increasing temperature. At large enough values of A, 
besides the maximum in x^^(0; T) there appears a mini- 
mum at a finite temperature which was also found for the 
isotropic Heisenberg chain [7] and contradicts the exact 
behavior [14]. Note that this artifact does not occur in 
the 2D XXZ model [13]. In the high-temperature limit, 
all susceptibilities depicted in Fig. 5 reveal a crossover to 
the Curie- Weiss behavior. 

Finally, in Fig. 6 the temperature dependence of the 
susceptibility Xn~ (^) ~ Xfi^ ("') is shown which again 
may be explained as SRO effect. Here, in the FM region, 
the transverse FM SRO results in a spin stiffness against 
the orientation of the transverse spin components along a 
staggered external field perpendicular to the ^-direction 
so that X^"!"") is suppressed at low temperatures and 
exhibits a max;imum. In the AFM region, the transverse 
AFM SRO results in an analogous temperature depen- 
dence of X7^~(0; T), where in the whole easy-plane region 
Tmax increases with increasing A. 

To summarize, we presented a Green's-function theory 
of magnetic SRO in the ID easy-plane XXZ model which 
allows, for the first time, the calculation of all static mag- 
netic properties over the whole easy-plane FM and AFM 
regions, whc^re no LRO occurs. That is, we computed 
the full wavenumber and temperature dependences of 
the anisotropic static spin susceptibilities and of the spin 
correlation functions of arbitrary range and at arbitrary 
temperature. In particular, in the FM region, we are able 
to reproduce the quantum-classical crossover in the lon- 
gitudinal spin correlations in qualitative agreement with 
the numerical diagonalization data of Ref. [2]. Moreover, 
the obtained maxima in the temperature dependences of 
the uniform longitudinal and transverse susceptibilities 
are explained as SRO effects. Prom the results of our 
theory we conclude that this approach may be STiccess- 
fuUy applied to other anisotropic spin models, where the 
description of spin correlations improves in higher dimen- 
sions [13]. 
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FIGURE CAPTIONS 

The authors would hke to thank Y. Gaididei, J. Stolze, 
and A. Weifie for helpful discussions. Fig. i. Transverse and longitudinal spin correlation 

functions (n = 1, 2) at T = 0. 



Fig. 2. Zero-temperature static susceptibilities at 
q = 0,w. 
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Fig. 3. Wavenumber dependence of the longitudi- 
nal (a) and transverse (b) static susceptibilities at T = 0. 

Fig. 4. Transverse correlation functions up to the 
fourth-nearest neighbors in the easy-plane ferromagnetic 
region. 

Fig. 5. Inverse uniform longitudinal (a) and staggered 
transverse susceptibility (b). 

Fig. 6. Temperature dependence of the uniform trans- 
verse susceptibility. 
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TABLE I. Longitudinal spin correlation functions C^^(r; A) for the easy-plane XXZ chain with A = ±0.3 at T = 0.1. The 
ED data given in parenthesis are taken from Ref . [2] . 



A 


n= 1 


n = 2 


n = 3 


n = 4 


0.3 


-1.15 X 10"°^ 


1.24 X 10"°^ 


-9.60 X 10"°^ 


1.34 X lO"""* 


-0.3 


-8.06 [-8.38] X 10"°^ 


-1.31 [-1.29] X 10"°^ 


-4.57 [-8.08] X 10""^ 


-1.89 [-2.70] X 10"°^ 



TABLE II. Temperature To (A; r) of the sign change in the longitudinal correlation functions C^^{T; A) of the ferromagnetic 
easy-plane XXZ model. The corresponding results obtained from ED of a 18-site ID chain [2] and a 2D 4x4 square lattice with 
periodic boundary conditions are given in parenthesis. 
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